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We discuss the salient features of the high density effective theory (HDET) of QCD, elaborating 
more on the matching for vector-vector correlators and axial-vector-vector correlators, which are 
related to screening mass and axial anomaly, respectively. We then apply HDET to discuss various 
color-superconducting phases of dense QCD. We also review a recent proposal to solve the sign prob- 
lem in dense fermionic matter, using the positivity property of HDET. Positivity of HDET allows 
us to establish rigorous inequalities in QCD at asymptotic density and to show vector symmetry 
except the fermion number is not spontaneously broken at asymptotic density. 



I. INTRODUCTION 



o 
o 

(N 

As physics advances, its frontier has expanded. One of the frontiers under active exploration is matter at extreme 
conditions. Recent surprising data, obtained from heavy-ion collisions and compact stars such as neutron stars, and 
also some theoretical breakthroughs have stimulated active investigation in this field . 

How does matter behave as we squeeze it extremely hard? This question is directly related to one of the fundamental 
questions in Nature; what are the fundamental building blocks of matter and how they interact. According to QCD, 
matter at high density is quark matter, since quarks interact weaker and weaker as they are put closer and closer. 

At what temperature and density does the phase transition to quark matter occur? To determine the phase diagram 
of thermodynamic QCD is an outstanding problem. The phases of matter are being mapped out by colliding heavy- 
jy-) , ions and by observing compact stars. Since QCD has only one intrinsic scale, Aqcd, the phase transition of QCD 
matter should occur at that scale as matter is heated up or squeezed down. Indeed, recent lattice QCD calculations 
found the phase transition does occur at temperature around 175 MeV 0. Even though lattice QCD has been 
quite successful at finite temperature but at zero density, it has not made much progress at finite density due to 
the notorious sign problem. The lattice calculation is usually done in Euclidean space and Euclidean QCD with a 
chemical potential has a complex measure, which precludes use of importance samplings, the main technique in the 
Monte Carlo simulation for lattice calculations. 
Oh, Lattice QCD at finite density is described by a partition function 

■ Z(p) =j dAdet{M)e~ s ^ A \ (1) 

where M = "f^Dg + jjr/g is the Dirac operator of Euclidean QCD with a chemical potential fi. The eigenvalues 
of M are in general complex, since is anti-Hermitian while [vy^ is Hermitian. For certain gauge fields such 

as A^(— x) = —A fl (x), M can be mapped into M' by a similarity transformation and thus its determinant M is 
nonnegative. However, for generic fields M ^ P~ X M^P and det (M) is complex. 

Recently there have been some progress in lattice simulation at small chemical potential, using a re-weighting 
method, to find the phase line Another interesting progress in lattice simulation was made at very high density 

in wnere it was shown that for QCD at high density the sign problem is either mild or absent, since the modes, 

responsible for the complexness of the Dirac determinant, decouple from dynamics or become irrelevant at high baryon 
density. 



II. HIGH DENSITY EFFECTIVE THEORY 



At low temperature or energy, most degrees of freedom of quark matter are irrelevant due to Pauli blocking. Only 
quasi-quarks near the Fermi surface are excited. Therefore, relevant modes for quark matter are quasi-quarks near 
the Fermi surface and the physical properties of quark matter like the symmetry of the ground state are determined 
by those modes. High density effective theory (HDET) 0,111 of QCD is an effective theory for such modes to describe 
the low-energy dynamics of quark matter. 

To find out the modes near the Fermi surface, one needs to know the energy spectrum of QCD, which is very 
difficult in general since it is equivalent to solving QCD. However, at high density /i » Aqcd, quarks near the Fermi 
surface carry large momenta and the typical interaction involves a large momentum transfer. Therefore, due to the 
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FIG. 1: Energy spectrum of quarks at high density 
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Fermi Surface 

FIG. 2: A patch covering the Fermi surface 



asymptotic freedom of QCD, the spectrum near the Fermi surface at high density looks very much like that of free 
fermion: (a. ■ p — /i + (3 m) i/)± = E±ip± , as shown in Fig. 1. We see that at low energy, E < 2/i, the states near the 
Fermi surface (\p\ ~ Pf), denoted as are easily excited while states deep in the Dirac sea, denoted as are 
hard to excite. 

At low energy, the typical momentum transfer by quarks near the Fermi surface is much smaller than the Fermi 
momentum. Therefore, similarly to the heavy quark effective theory, we may decompose the momentum of quarks 
near the Fermi surface as 



p» = /i^ + Z", 



(2) 

where = (0,vf) and vf is the Fermi velocity. For quark matter, the typical size of the residual momentum is 
~ Aqcd, and the Fermi velocity of the quarks does not change for /x 3> Aqcd, when they are scattered off by soft 
gluons. 

We now introduce patches to cover the Fermi surface, as shown in Fig. 2. The sizes of each patch are 2A in vertical 
direction to the Fermi surface and 2Aj_ in horizontal direction. The quarks in a patch are treated to carry a same 
Fermi velocity. 

The energy of the quarks in the patch is given as 



E = -n+ yjp 1 + m 2 = I ■ v F + — + O 



2/i 



(3) 



We see that at the leading order inl/fi expansion, the energy is independent of the residual momentum, l±, perpen- 
dicular to the Fermi velocity. In HDET, therefore, the perpendicular momentum labels the degeneracy and should 
satisfy a normalization condition 



J2 / d 2 f_L = 47Tp£. 



patches 1 

To identify the modes near the Fermi surface, we expand the quark field as 
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FIG. 3: Tree-level matching: The double line denotes ip- modes and the single line tp + . 



where ip±(vp,x) satisfies respectively 

1 ± a ■ vf 



ijj± = %/j±. (6) 



Note that the projection operator P± = (1 ± a. ■ vf)/2 projects out the particle state, ip+, and the anti-particle state, 
■0_ (or more precisely ip-), from the Dirac spinor field ^ . The quasi-quarks in a patch carries the residual momentum 
and is given as 

i>+(v F ,x)= l + d 2 VF e'^^(x) (7) 
The Lagrangian for quark fields becomes 

£ = § (ij> + ^7°) * = £ $ (p + + P _) ( fl y +p) ( P+ + p_)^ 

VF 

= V+ + V>- ( 2a*7° + + VP- if>± -P+ + h.c] , (8) 

where we neglected the quark mass term for simplicity and = (1,vf)- The parallel component of the covariant 
derivative is D 1 ^ = D ■ V and the perpendicular component D± = D — D\\. From the quark Lagrangian one can 
read off the propagators for ^±(^,2): 



2/i 



(9) 



We see indeed that in HDET the quarks near the Fermi surface or ip + modes are the propagating modes, while ip_ 
are not. 

By integrating out ?/>_ modes and the hard gluons, one obtains the high density effective theory of QCD. In general 
the integration results in nonlocal terms in the effective theory and one needs to expand them in powers of \j [i. This 
is usually done by matching the one-light-particle irreducible amplitudes of the microscopic theory with those of the 
effective theory. For tree-level amplitudes, this is tantamount to eliminating the irrelevant modes, using the equations 
of motion. 

For instance, a one-light particle irreducible amplitude in QCD of two gluons and two quarks is matched as 

P ± iP-(v F ,x)$-i P ± -tP+(v F ,y) = $+iP± (-^7) * P±*l>+ - ( U ) 

which is shown in Fig. 3. Similarly one can eliminate the hard gluons. Integrating out hard gluons results in four-Fermi 
interactions of tp+ modes. (See Fig. 4.) One continues matching one-loop or higher- loop amplitudes. One interesting 
feature of HDET is that a new marginal operator arises at the one-loop matching, when incoming quarks are in 
Cooper-paring kinematics, namely when they have opposite Fermi velocities. As shown in Fig. 5, when the incoming 
quarks have opposite Fermi velocity, the amplitudes in HDET are ultra-violet divergent while QCD amplitudes are 
not. Therefore, one needs to introduce a four-Fermi operator as a counter term to remove the UV divergence. If we 
collect all the terms in the effective theory, it has a systematic expansion in and coupling constants a s as 

£HDET = 6 1 ^ + z 7 p^ + -|i^ +7 (^) 2 ^ + + (12) 

where b\ = 1 + 0(a s ), ci = 1 + 0(a s ), ■ ■ ■. Note that HDET has a reparametrization invariance, similarly to 
heavy quark effective theory, which is due to the fact that the Fermi velocity of quarks in a patch is not uniquely 
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determined. For a given quark momentum, the corresponding Fermi velocity is determined up to reparametrization; 
Hp — > vp+Sl±/ /j, and I — > I— SI, where 61 ± is a residual momentum perpendicular to the Fermi velocity. As in the heavy 
quark effective theory Q , the renormalization of higher-order operators are restricted due to the reparametrization 
invariance. For instance, b\ — c\ at all orders in a s . 

In order for the effective theory to be meaningful, it should have a consistent power-counting. We find the consistent 
counting in HDET to be for = A 

To be consistent with the power counting, we impose in loop integration 

12 



d 2 l± 11 = for n > 0. (14) 

So far we have restricted ourselves to operators containing quarks in the same patch. For operators with quarks in 
different patches, one has to be careful, since the loop integration might jeopardize the power-counting rules. Indeed, 
consistent counting is to sum up all the hard- loops, as shown by Schafer |lCj . 



III. MORE ON MATCHING 



In HDET, the currents are given in terms of particles and holes but without antiparticles as 

J» = J2$(VF,x)^($F,x)-^[rt,ip X }iP + --- , (15) 

where the color indices are suppressed and we have reverted the notation ip for ip + henceforth. We find that the 
HDET current is not conserved unless one adds a counter term. Consider the current correlator 

(J»(x)r(y)) = - f T f = [ e Hp-(x-v)nM"(p) (16) 
6A tl [x)6A v {y) J p 

where the vacuum polarization tensor 

nTM = -^f d ^( -%f V : V : ) (17) 

2 J 47r \p ■ V + iep ■ vp I 
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FIG. 6: Matching two-point functions 



and M 2 — N f g 2 fi 2 / (2tt 2 ) . We see that the vacuum polarization tensor is not transverse, p^Yl^ijp) 7^ 0, which means 
that the current is not conserved. The physical reason for this is that not only modes near the Fermi surface but 
also modes deep in the Fermi sea respond to external sources collectively. To recover the current conservation in the 
effective theory, we need to add the DeBye screening mass term due to ip_ (See Fig. 6): 



-pcff 1 ^ -pcff j^eff 



M 2 



Then vacuum polarization tensor becomes 



= 11"" 



Vp 



(18) 



(19) 



The modified polarization tensor is now transverse, p^IT 1 " = 0. 

Now, let us consider the divergence of axial currents in HDET, which is related to the axial anomaly and also to 
how the quark matter responds to external axial-current sources like electroweak probes. 

It is easy to show that the axial anomaly in dense matter is independent of density or the chemical potential \x |ll| . 
In general one may re- write the divergence of axial currents in dense QCD as follows: 



16tt 



(20) 



where the first term is the usual axial anomaly in vacuum and the second term is due to matter. However, one can 
explicitly calculate the second term, which is finite, to find A a/3 (fj,) = 0. In HDET, the axial anomaly due to modes 
near the Fermi surface is given as 



E / e ik ^ +ik ^ (d^(v F ,0)J a (v F ,x)jP(v F ,y)) = A^ 

- J x.y 

By explicit calculation we find 

= • \ (£i x k 2 ) 1 , A% = (k 10 k 2l - k u k 20 ) ■ 



(21) 



(22) 



We see that the modes near the Fermi surface contributes only some parts of the axial anomaly. As in the vector 
current, the rest should come from modes in the deep Fermi sea and from anti-particles. To recover the full axial 
anomaly we add a counter term (See Fig. 7), which is two thirds of the axial anomaly plus a Chern-Simons term: 



A o ^ = + £^e a(3p °k lp k 2a + ^e a(301 (k w k 2l - k u k 20 ) . 



(23) 



IV. COLOR SUPERCONDUCTIVITY IN DENSE QCD 



At high density, quarks in dense matter interact weakly with each other and form a Fermi sea, due to asymptotic 
freedom. When the energy is much less than the quark chemical potential {E <C /i), only the quarks near the Fermi 
surface are relevant. The dynamics of quarks near the Fermi surface is effectively one-dimensional, since excitations 
along the Fermi surface do not cost any energy. The momentum perpendicular to the Fermi momentum just labels 
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FIG. 7: Matching axial anomaly. + denotes ip+ and — denotes ip-. 



the degeneracy, similarly to the perpendicular momentum of charged particle under external magnetic field. This 
dimensional reduction due to the presence of Fermi surface makes possible for quarks to form a Cooper pair for any 
arbitrary weak attraction, since the critical coupling for the condensation in (1+1) dimensions is zero, known as the 
Cooper theorem in condensed matter. 

While, in the BCS theory, such attractive force for electron Cooper pair is provided by phonons, for dense quark 
matter, where phonons are absent, the gluon exchange interaction provides the attraction, as one-gluon exchange 
interaction is attractive in the color anti-triplet channel [3ll| One therefore expects that color anti-triplet Cooper pairs 
will form and quark matter is color superconducting, which is indeed shown more than 20 years ago |l3l Il4l Il5j| . 

At intermediate density, quarks and gluons are strongly interacting and gluons are therefore presumably screened. 
Then, QCD at intermediate density may be modelled by four-Fermi interactions and higher-order terms by massive 
gluons. 
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£qcd 9 -^H>H> + • • • , (24) 

where the ellipsis denotes higher-order terms induced by massive gluons. When the incoming quarks have opposite 
momenta, the four-Fermi interaction is marginally relevant, if attractive, and all others are irrelevant. As the renor- 
malization group flows toward the Fermi surface, the attractive four-Fermi interaction is dominant and blows up, 
resulting in a Landau pole, which can be avoided only when a gap opens at the Fermi surface. This is precisely 
the Cooper-instability of the Fermi surface. The size of gap can be calculated by solving the gap equation, which is 
derived by the variational principle that the gap minimizes the vacuum energy: 

AW A) A r d±k A 

8 A G J {2^Ykl-{k-v F f-^ 1 1 



which gives 



d 4 k 

W [(1 + it)k Q f - (k ■ v F Y 



We note that the integrand in Eq. does not depend on k± , whose integration gives the density of states at the 
Fermi surface, and the it prescription is consistent with the Feynman propagator. The pole occurs at 



k Q = ±\J {k ■ v F f + A 2 Tie (27) 
or in terms of full momentum p = fx v + k it occurs at 

p Q = ±v/(b1-M) 2 + A 2 T *e • (28) 

We find the solution to the gap equation 

A„ = 2Ai exp(-^). (29) 

For generic parameters of dense QCD, the gap is estimated to be 10 ~ 100 MeV at the intermediate density. The free 
energy of the BCS state is given as 

= A ~§- f\x + 9 2 \nx)d X = -^- 2 Al, (30) 




FIG. 8: Eliashberg equation at high density. 

where x = A /(2m) and g 2 = 2G/j, 2 /ir 2 . At high density magnetic gluons are not screened though electric gluons are 
screened 0, 03, H] • The long-range pairing force mediated by magnetic gluons leads to the Eliashberg gap equation 
(See Fig. 8). 

Abo = / dq =ln ] r , (31) 

where A = • (^/M) 5 e 3 / 2 « and £ is a gauge parameter. Due to the unscreened but Landau-damped gluons, there 
is an extra (infrared) logarithmic divergence in the gap equation, when the incoming quark momentum is collinear 
with the gluon momentum. The Cooper-pair gap at high density is found to be 0, QjJ 
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The numerical prefactor of the gap is not complete, since the contributions from subleading corrections to the gap 
equation that lead to logarithmic divergences, such as the wavefunction renormalization and the vertex corrections, 
are not taken into account. Recently, however, the contributions to the prefactor, coming from the vertex corrections 
and the wavefunction renormalization for quarks were calculated by finding a (nonlocal) gauge poj . where the quark 
wavefunction is not renormalized for all momenta, Z(p) = 1. At the nonlocal gauge, £ ~ 1/3. The subleading 
corrections therefore increase the leading-order gap at the Coulomb gauge by about two thirds. 



V. QUARK MATTER UNDER STRESS 



It is quite likely to find dense quark matter inside compact stars like neutron stars. However, when we study the 
quark matter in compact stars, we need to take into account not only the charge and color neutrality of compact stars 
and but also the mass of the strange quark, which is not negligible at the intermediate density. By the neutrality 
condition and the strange quark mass, the quarks with different quantum numbers in general have different chemical 
potentials and different Fermi momenta. When the difference in the chemical potential becomes too large the Cooper- 
pairs breaks or other exotic phases like kaon condensation or crystalline phase is more preferred to the BCS phase. 

Let us consider for example the pairing between up and strange quarks in chemical equilibrium. The energy 
spectrum of up quarks is given as 

E = -fi±\p\, (33) 

while the energy of strange quarks of mass M s becomes 

E = -^±^\p\ 2 + M 2 . (34) 

The Fermi sea of up and strange quarks is shown in Fig. 9. Because of the strange quarks mass, they have different 
Fermi momenta. Note that the Cooper-pairing occurs for quarks with same but opposite momenta. Therefore, at 
least one of the pairing quarks should be excited away from the Fermi surface, costing some energy. Let us suppose 
that the Cooper-pair gap opens at \p\ = p between two Fermi surfaces, p s F < p < 

To describe such pairing, we consider small fluctuations of up and strange quarks near p. The energy of such 
fluctuations of up and down quarks is respectively 

E u = -n+\P + T\ - Sfi u + v u -l, E s ~ -Sn s +v s -l, (35) 

where Sfj, u = ^ — p and 5[i s = [X — \Jp 2 + M 2 . v u and v s are the velocities of up and strange quarks at \p] = p. Let 
A be the BCS gap for the u, s pairing. Then, the Lagrangian for the u, s quarks is given as 

C = u{i ft + /i7°) u + s c (i ft — /i7° — M s ) s c — As c u + h.c. + d nt , (36) 
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FIG. 9: Fermi sea of up and strange quarks. 



where s c is the charge conjugate field of s quark. In HDET, the Lagrangian becomes 

Ajdet 3 (iV u ■ d + 8ll u ) u + sl (iV s ■ d - Sll s ) s c - As c u + h.c, 
where V u = (l,v u ) and V s — (1, — v s ). The Cooper-pair gap equation is then 

iA(l)K(p-l) 



A(p) = 



(1 +ie)l - l-v u + 8li u (1 +ie)l + I ■ v s -6/j, u 



-A 2 



(37) 



(38) 



where K is the kernel for the gap equation and is a constant for the four-Fermi interaction. By examining the pole 
structure, we see that the Cooper-pair gap does not exist when 



A 2 



(39) 



Only when — 5fi u Sfi s < A 2 /4, one can shift Iq — > 1' = lo + Sfi u or Iq — + 1' = Iq — <5/i s without altering the pole structure. 
Note that the gap becomes biggest when Sfi u = —5[i s (= 8 /a), which determines the pairing momentum to be 



p = fj, 



Ml 
4/z ' 



(40) 



If 8li < A/2 or A > Af 2 /(2/i), the solution to the Cooper-pair gap exists. The gap equation then can be written as, 
shifting Iq, in Euclidean space 



A(p) 



dH A(0 



(2^) 4 q + a 2 



K(l-p), 



(41) 



where I?, = Iq + c 2 (l ■ v) 2 and c 2 = p/y^p 2 + M 2 . In HDET, one can easily see that the Cooper-pair gap closes if 
the effective chemical potential difference, 2<5/i, due to an external stress, exceeds the Cooper-pair gap when there 
is no stress. One should note that even before the Cooper-pair gap closes other gap may open as shown by many 
authors 0,H2. But, one needs to compare the free energy of each phases to find the true ground state for quark 
matter under stress. 



VI. POSITIVITY OF HDET 



Fermionic dense matter generically suffers from the sign problem, which has thus far precluded lattice simula- 
tions |23l ]. However, the sign problem usually associated with fermions is absent if one considers only low-energy 
degrees of freedom. The complexness of the measure of fermionic dense matter can be ascribed to modes far from 
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the Fermi surface, which are irrelevant to dynamics at sufficiently high density in most cases, including quark mat- 
ter [HEl- For modes near the Fermi surface, there is a discrete symmetry, relating particles and holes, which pairs the 
eigenvalues of the Dirac operator to make its determinant real and nonnegative. Especially, the low energy effective 
theory of dense QCD has positive Euclidean path integral measure, which allows one to establish rigorous inequalities 
that the color-flavor locked (CFL) phase is the true vacuum of three flavor, massless QCD. 

As simple example, let us consider a fermionic matter in f+f dimensions, where non-relativistic fermions are 
interacting with a gauge field A. The action is in general given as 

drdx ift [{-d T + # + e F ) - e(-id x + A)) </>, (42) 



where e(p) ~ p 2 /(2m) + • • ■ is the energy as a function of momentum. Low energy modes have momentum near the 
Fermi points and have energy, measured from the Fermi points, 



dE 

E(p±p F ) ~±v F p, v F = — 

op 



(43) 



If the gauge fields have small amplitude and are slowly varying relative to scale p F , the fast modes are decoupled 
from low energy physics. The low energy effective theory involving quasi particles and gauge fields has a positive, 
scmi-dcfinitc determinant. 

To construct the low energy effective theory of the fermionic system, we rewrite the fermion fields as 

1>(x,t) = i> L (x,T)e +ipFX + Mxrfe-*"', (44) 

where ipL.R describes the small fluctuations of quasiparticles near the Fermi points. Using e ±WFX E(—id x + 
A) e TipFX \p{x) « ± v F {-id x + A)ip(x), we obtain 

Scft= U\{-d T + i<f) + id x -A)^ L +^{-d T + i<j)~id x +A)^ R \. (45) 

Jt,x 

Introducing the Euclidean (1+1) gamma matrices 70.1,2 and ipL.R = 5(1 ± ^2)^, we obtain a positive action: 

S eS = I drdx^idf.+iA^j^ [ drdx^p^. (46) 



Since Jp= 72^72, the determinant of _^is positive, semi-definite. 

In this example, we see that near the Fermi surface, modes have low energy, slowly varying, and thus lead to an 
effective theory without any sign problem what so ever if they couple to slowly varying background fields. QCD at 
high baryon density falls into this category, since the coupling constant is small at high energy due to asymptotic 
freedom. 

HEDT of quark matter is described by 

£hdet - <M7p^+ ~ ^+7°U\) 2 V>+ + (47) 
where 7^ = (7 , v F v F ■ 7) = 7^ — 7^. We see that the leading term has a positive determinant, since 

M eft = lf-D(A) = 75^75. (48) 

In order to implement this HDET on lattice, it is convenient to introduce an operator formalism, where the velocity 
is realized as an operator, 

"= i=i ■ (49) 

Then the quasi quarks near the Fermi surface become 

rp — exp (+i[ix ■ v) — ip+. (50) 

Now, neglecting the higher order terms, the Lagrangian becomes with X — exp(ifix ■ v)(l + a • v)/2, 

£ H DET=^+7f (^+a+)</>+, (5!) 
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FIG. 10: Simulation with two lattices with different lattice spacings 



where = A M X denotes soft gluons whose momentum < fi. Since v ■ dv • 7 = d • 7 , we get 



(52) 



which shows that the operator formalism automatically covers modes near the full Fermi surface. 

Integrating out the fast modes, modes far from the Fermi surface and hard gluons, the QCD partition function Q 
becomes 



Z(n) = / d,4 + det (Mea) e 



(53) 



where 



_pa pa 

4 \iv 



M 2 

AO. A a 



(54) 



and A± = A — Ah, the Debye mass M = \J Nf / (2n 2 )g s fi . At high density the higher order terms (~ A//x) are 
negligible and the effective action becomes positive, semi-definite. Therefore, though it has non-local operators, 
HDET in the operator formalism, free from the sign problem, can be used to simulate the Fermi surface physics 
like superconductivity. Furthermore, being exactly positive at asymptotic density, HDET allows to establish rigorous 
inequalities relating bound state masses and forbidding the breaking of vector symmetries, except baryon number, in 
dense QCD [|. 

With the help of previous two examples, we propose a new way of simulating dense QCD, which evades the sign 
problem. Integrating out quarks far from the Fermi surface, which are suppressed by 1/fi at high density, we can 
expand the determinant of Dirac operator at finite density, 



det (M) = [real, positive] 



1 



(55) 



As long as the gauge fields are slowly varying, compared to the chemical potential /1, the sign problem can be evaded. 
As a solution to the sign problem, we propose to use two lattices with different spacings, a finer lattice with a lattice 
spacing a^ct ~ (J> f° r fcrmions and a coarser lattice with a lattice spacing a gaug0 <C /i -1 for gauge fields and then 
compute the determinant on such lattices. 

The determinant is a function of plaquettes {U XM } which are obtained by interpolation from the plaquettes on the 
coarser lattice of lattice spacing a gaugo . To get the link variables for the finer lattice, we interpolate the link variables 
U XM G SU(3) (see Fig. 1100 : Connect any two points 171,32 on the group manifold as 



g(t) = gi + t(g 2 -gi) , 0<t<l. 



(56) 



For importance sampling in the lattice simulation, one can use the leading part of the determinant, [real, positive]. 
This proposal provides a nontrivial check on analytic results at asymptotic density and can be used to extrapolate 
to intermediate density. Furthermore, it can be applied to condensed matter systems like High-T c superconductors, 
which in general suffers from a sign problem. 

Positivity of the measure allows for rigorous QCD inequalities at asymptotic density. For example, inequalities 
among masses of bound states can be obtained using bounds on bare quasiparticlc propagators. One subtlety that 
arises is that a quark mass term does not lead to a quasiparticle gap (the mass term just shifts the Fermi surface). 
Hence, for technical reasons the proof of non-breaking of vector symmetries [24| must be modified. (Naive application 
of the Vafa-Witten theorem would preclude the breaking of baryon number that is observed in the color-flavor-locked 
(CFL) phase [25|l. A quasiparticle gap can be inserted by hand to regulate the bare propagator, but it will explicitly 
violate baryon number. However, following the logic of the Vafa-Witten proof, any symmetries which are preserved 
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by the regulator gap cannot be broken spontaneously. One can, for example, still conclude that isospin symmetry is 
never spontaneously broken (although see below for a related subtlety) . In the case of three flavors, one can introduce 
a regulator d with the color and flavor structure of the CFL gap to show rigorously that none of the symmetries of 
the CFL phase are broken at asymptotic density. On the other hand, by applying anomaly matching conditions |26| . 
we can prove that the SU(3)a symmetries are broken. We therefore conclude that the CFL phase is the true ground 
state for three light flavors at asymptotic density, a result that was first established by explicit calculation 0, l27t l28j . 

To examine the long-distance behavior of the vector current, we note that the correlator of the vector current for 
a given gauge field A can be written as 



(j;(v F ,x)J b u (v F ,y)) = -Tr ^T a S A (x,y;dh„T b S A (y,x;d), 



(57) 



where the SU(Nf) flavor current J^(vf,x) — ip+(i)F, x)^^T a if}+{vF ^ x). The propagator with S'f7(3)v r -invariant IR 
regulator d is given as 



1 r°° 

S A (x, y; d) = (x\ Jj \y)= dr (x\ e -"C-**0 {y) 



where with D = d + iA 



M = 70 



D ■ V d 
dt D ■ V 



Since the eigenvalues of M are bounded from below by d, we have 

1 



(*l m M 



,-dR 



Jr 



■y/(x\x)y/(y\y), 



(58) 



(59) 



where R = \x — y\. The current correlators fall off rapidly as R — > oo; 



d^et M eff (A) e" s °« { J A (v F , x) J? (v F , y) ) 



< 



(J, 



A/„-t 



X)J„(VF; 



y)Y 



< 



-2dR 



f/ 2 



\(x\x)\\(y\y)\, 



(60) 



where we used the Schwartz inequality in the first inequality, since the measure of the effective theory is now positive, 
and equation l|59|) in the second inequality. The IR regulated vector currents do not create massless modes out of 
the vacuum or Fermi sea, which implies that there is no Nambu-Goldstone mode in the SU(3)v channel. Therefore, 
for three massless flavors SU(3)v has to be unbroken as in CFL. The rigorous result provides a non-trivial check on 
explicit calculations, and applies to any system in which the quasiparticle dynamics have positive measure. 

It is important to note the order of limits necessary to obtain the above results. Because there are higher-order 
corrections to the HDET, suppressed by powers of A//i, that spoil its positivity, there may be contributions on the 
RHS of of the form 



0[±) /(./?.. 



(61) 



where f(R) falls off more slowly than the exponential in l|60|l . To obtain the desired result, we must first take the 
limit /i — ► oo at fixed A before taking R — * oo. Therefore, our results only apply in the limit of asymptotic density. 

Although our result precludes breaking of vector symmetries at asymptotic density in the case of three exactly 
massless quarks [3(3, it does not necessarily apply to the case when the quark masses are allowed to be slightly 
non-zero. In that case the results depend on precisely how the limits of zero quark masses and asymptotic density 
are taken, as we discuss below. 

In the authors investigate the effect of quark masses on the CFL phase. These calculations are done in the 
asymptotic limit, and are reliable for sufficiently small quark masses. When m u = rrid = m << m s (unbroken SU(2) 
isospin, but explicitly broken SU(3)), one finds a kaon condensate. The critical value of m s at which the condensate 
forms is m* ~ m 1 / 3 A^ 3 , where Ao is the CFL gap (see, in particular, equation (8) of the first paper). As kaons 
transform as a doublet under isospin, the vector SU(2) symmetry is broken in seeming contradiction with our result. 

However, a subtle order of limits is at work here. For simplicity, let us set m = 0. Note that the CFL regulator 
d, which was inserted by hand, explicitly breaks SU(3)a through color-flavor locking, leading to small positive mass 
squared for the pions and kaons, given as 



a s d In 



(S) 



(62) 
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The meson mass is not suppressed by 1//J., since, unlike the Dirac mass term, the regulator, being a Majorana mass, 
does not involve antiquarks [23 . 

Therefore, even when the light quarks are massless, there is a critical value of m s necessary to drive negative the 
mass-squared of kaons and cause condensation: 

1/2 



g^/xlng)] > (dn) 1/2 , (63) 



where g s is the strong coupling constant. Note the product of g s with the logarithm grows as /i gets large. To obtain 
our inequality we must keep the regulator d non-zero until the end of the calculation in order to see the exponential 
fall off. To find the phase with kaon condensation identified in we must keep m s larger than m*. (Note fx — > 00, 
so to have any chance of finding this phase we must take d — > keeping dR large and d[i small.) 
Since the UV cutoff of the HDET must be larger than m s , we have 



which implies 



i> if) >|£. m 



±f(R) > jf(R). (65) 



Note the right hand side of this inequality does not necessarily fall off at large R, and also does not go to zero for 
fx — > 00 at fixed A and d. This is a problem since to apply our inequality the exponential falloff from l|60|l must 
dominate the correction term (|61[) . which is just the left hand side of (|65|l . Combining these equations, we see that 
the exponential falloff of the correlator is bounded below, 

-2dR j 

> t/(«)> (66) 



d 2 A 

in the scaling region with a kaon condensate, m s > m*. 

Alternatively, if we had taken m s to be finite for fixed regulator d (so that, as /1 — + 00, eventually m s < m*), the 
inequality in ijoTl could be applied to exclude a Nambu-Goldstone boson, but we would find ourselves in the phase 
without a kaon condensate. 
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